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Disclaimer 


•  You  will  not  receive  any  personal  financial 
advice  during  this  talk,  as  I  am  not 
officially  qualified  or  certified  to  do  so. 

•  However,  my  presentation  is  intended  to 
get  you  to  think  mathematically  about  one 
of  the  retirement  savings  options  available 
to  many  of  you. 
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Questions  for  Consideration 


•  Why  might  I  be  more  risk  tolerant  than  I 
currently  believe? 

•  What  are  the  L  (Lifecycle)  funds?  How  are 
they  constructed?  Why  might  they  be  of 
interest  (or  not)  to  me? 

•  What  if  stock  and  index  fund  returns  are  not 
normally  distributed,  as  is  commonly 
assumed? 

•  How  does  the  choice  of  reward  and  risk 
measures  affect  optimal  TSP  portfolios? 
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Big  Picture  on  Saving  $$$ 


“There  is  no  scholarship  for  retirement!” 


-Unknown 

Spouse’s  401  (k)  with  matching  funds 
Roth  IRA  (for  Soldier/civilian  and  spouse) 
Thrift  Savings  Plan  (TSP) 

Spouse’s  401  (k)  without  matching  funds 
Coverdale  Educational  Savings  Accounts 
529  Tuition  Plans  (prepaid  or  savings) 

(ordering  of  these  depends  on  tax  considerations) 


tTy. 


Jr 


Robert  H.  Smith 

£  SCHOOL  OF  BUSINESS 


©  2008  Scott  T.  Nestler 


Leaders  for  the  Digital  Economy 


Thrift  Savings  Plan  (TSP)  Overview 


•  Largest  defined  contribution  retirement 
savings  and  investment  plan 

-3.7  million  participants 
-$210  billion  in  assets 

•  401  (k)  equivalent  for  government 
employees  and  uniformed  service 
members 

•  5  non-traded  core  funds 

•  Can  rebalance  daily  with  no  direct  costs 
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RISKY”  ASSETS 


Core  TSP  Funds 


Fund 

Description 

Assets* 

Mean 

Return# 

Standard 

Deviation# 

G 

short-term,  specially  issued 

Treasury  securities 

$66. 6B 
(39.2%) 

6.4% 

1.5% 

F 

tracks  Lehman  Brothers  U.S. 
Aggregate  (LBA)  Index 

$10. 2B 
(6.0%) 

7.3% 

5.6% 

c 

tracks  S&P  500  Index 

$66. 7B 
(39.3%) 

13.0% 

17.9% 

s 

tracks  Dow  Jones  Wilshire  4500 
Completion  Index 

$13. 7B 
(8.1%) 

13.3% 

19.9% 

1 

V _ * 

tracks  MSCI  EAFE  (Europe, 
Australia,  Far  East)  Index 

$12. 6B 
(7.4%) 

7.8% 

18.7% 

*  As  of  Dec  31 , 2005  *  For  the  period  1 988-2005 
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Returns  /  Investment  Horizon 


•  Returns 

-  Arithmetic: 
-Log  : 


r. 

i,t 


c  _  e 

^i,t 


r. 

i,t 


=  ln(Sw)-ln(S/jM)  =  ln 


S,,  ' 


Su- . 


y 


•  Investment  Horizon:  20  years 

-  Point  of  military  (not  ultimate)  retirement 

-  System  encourages  20  year  careers 

-  Employment  options  vary  greatly 

-  Can  move  TSP  assets  to  other  plans 
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Daily  Returns  Time  Series 
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2000  4000  6000  8000  10000  12000 


Daily  Returns  Distributions 


G  Fund 


0.00010  0.00015  0.00020  0.00025 


C  Fund 


•  G  fund  appears  approximately  Gaussian 

•  C,  S,  I,  and  F  funds  are  more  peaked  with  heavy  tails 

•  Goodness  of  Fit  testing  at  common  levels  of  significance 
rejects  Normal  for  F,  C,  and  S  funds,  even  with  batched  means 

Mr  Robert  H. Smith 
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Mean-Variance  Portfolio  Optimization 

(Markowitz,  1952) 


Let  Xj  =  fraction  of  funds  invested  in  asset  i 
Rt  =  expected  return  of  asset  i 
Rp  =  expected  return  of  portfolio  p 
of  =  variance  of  return  of  asset  i 
<jjk  =  covariance  of  return  of  asset  j  with  asset  k 


N 


N  N 


Minimize  a2P=^  X)cj)  +LZ  XjX^jk 


7=1  k= 1 
k*j 


N 


Subject  to:  Zv=i 


/=i 

N 
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L  (Lifecycle)  Funds 

Invest  in  5  core  TSP  funds  based  on  time  horizon  to 
provide  highest  possjble  rate  of  return  for  risk  taken. 


C 

O 

+■» 

TO 

U 

o 


30  28  26  24  22 


L2040 


L2030  L2020 

Years  Until  Retirement 


•  Over  time,  investments  shift  away  from  stocks  and 
into  bonds. 

•  L  Funds  are  great,  but ... 
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Reward-Risk  Profile  of  TSP  Funds 
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VG-ICA  Factor  Model 

(Madan  &  Yen ,  2004) 

(R  —  n )  —  XB  +  g 
D  =  XB  +  g 

1 .  Use  Independent  Component  Analysis 
(ICA)  on  asset  returns  D  to  identify 
underlying  factors  X 

2.  Fit  the  Variance  Gamma  (VG)  distribution  to 
each  retained  factor  by  MLE 

3.  Use  Expected  Utility  to  determine  optimal 
portfolio  of  VG-ICA  factors;  convert  back  to 
optimal  portfolio  of  assets  (TSP  funds) 

Mr  Robert  H. Smith 
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Independent  Component  Analysis  (ICA) 


(Underlying  (Daily  log  (Estimates  of 


Factors)  Returns)  Factors) 

•  Principal  Component  Analysis  (PCA) 

-  Focus  on  finding  uncorrelated  components  in  Gaussian  data 

-  Maximizes  explained  variance 

-  Uses  second-order  statistics 

•  Factor  Analysis 

-  Essentially  PCA  with  extra  terms  to  model  noise 

•  ICA 

-  Focus  on  independent  and  non-Gaussian  components 

-  Maximizes  non-Gaussianity  (to  maximize  information) 

-  Uses  higher-order  statistics 

Mr  Robert  H. Smith 
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-1.0  0.0  1.0  -1.0 


Another  ICA  Example 

Original  Signals  (s)  Mixed  Signals  (x)  ICA  source  estimates  (y) 
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ICA  versus  PCA 


■  Principal  Component  Analysis  (PCA)  finds: 

■  directions  of  maximal  variance  in 
Gaussian  data  (second-order  statistics). 


■  directions  of  maximal  variance  in  non- 
Gaussian  data  (second-order  statistics). 


■  Independent  Component  Analysis  (ICA)  finds 
directions  of  maximal  independence  in  non- 
Gaussian  data  (higher-order  statistics). 
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How  Many  ICs  to  Keep? 

“Scree”  Plot  R2  Values  from  Regression 


#  ICs  F 

Kept  Fu 

:  C 

nd  Fund 

S 

Fund 

1 

Fund 

5 

1  1 

1 

1 

4  ,9S 

180  .9994 

.9998 

.9949 

3 

.8983 

.9652 

.9869 

2 

.8962 

.9629 

.127# 

1  2  3  4  5 

ICA  Factor  Number 

•  Dropping  more  than  one  1C  reduces  fit  on  at  least 
one  fund 

•  The  first  four  have  excess  kurtosis 
— »  Keep  4  Independent  Components  (ICs) 
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VG  Process  and  Distribution 


Pure  jump  process  with  two  representations 

-  Time-changed  Brownian  motion  (Madan  &  Seneta,  1990) 

XVG  (t;  v,  0,  cr)  =  b(y(t;  1,  v),  0,  <j) 

-  Difference  of  2  Gamma  processes  (Madan,  Carr  &  Chang,  1998) 

XrG(t)  =  Gp(t)-Gn(t) 

Parameters:  a  controls  spread 

v  affects  kurtosis 
0  impacts  skewness 


Density  Function  (Madan,  Carr  &  Chang,  1998) 


h{z)  = 


2  exp  (Ox  /  <t2  ) 


f 


v'lv  42n<jT{-) 

V 


X 


j__y 
\lv  4 


V 


2  a2/v  +  0z 


J 


K{  ,  /v  +  02) 

v~2 


with  x  =  z-6  where  z  =  ln(S ( t ) / S(t- 1)) 
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Examples  of  VG  Distributions 


Effect  of  v  Effect  of  0 


x 
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Fitting  VG  by  MLE 

•  Given  observed  IID  data  Xp  X2,...,  Xn,  define  the 

likelihood  function  as: 

•  The  MLE  (maximum  likelihood  estimator)  0 

maximizes  L(0)  over  all  permissible  values  of  0. 

•  Actually,  maximizing  the  log  likelihood  function 

/n(L(0))  is  easier 

•  For  the  VG  distribution  with  three  parameters,  this 
becomes: 

K<T,V,0)  =  In  L(c 7,  v,  (?)  =  '£  (X,) 

i 

(using  pdf  given  before  from  Madan,  Carr,  and  Chang,  1998) 
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Comparison  of  Fitted  VG  and  Normal(0,1) 

IC1  IC2 


FactyCQeturn 


-4  -2  0  2  4 

Faclj^jjteturn 


r 

<? 


Factor  Return 

Robert  H.  Smith 

SCHOOL  OF  BUSINESS 


Factor  Return 


o  Data 

- Normal(0,1) 

.  Fitted  VG 
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Fitted  VG  Parameters  /  Chi-Square  Statistics 


IC# 

Fitted  VG  Parameters  -  Daily 
(Annualized) 

X2  Test  Statistic  (p-values) 
(X2,0i,17=  33.41) 

<7 

V 

e 

VG(a,v,0) 

N(0,1) 

IC1 

0.933 

(14.814) 

0.969 

(.00385) 

-0.015 

(-3.773) 

IIIIIIIIH1IIIIII 

546.89 

(7.83E-104) 

IC2 

0.980 

(15.558) 

0.820 

(.00326) 

-0.001 

(-0.222) 

13.37 

(0.82) 

295.73 

(1.50E-51) 

IC3 

0.989 

(15.703) 

0.586 

(.00232) 

-0.004 

(-1.019) 

34.94 

(0.014) 

171.70 

(1.32E-26) 

IC4 

0.991 

(15.739) 

0.468 

(.00186) 

-0.005 

(-1.149) 

25.29 

(0.151) 

118.07 

(2.54E-16) 

•  Some  excess  kurtosis  and  slight  negative 
skewness  in  each  1C 

•  VG  fits  much  better  than  Normal  distribution 

Mr  Robert  H. Smith 
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Utility  Theory  &  Risk  Aversion 

Utility-  a  measure  of  relative  satisfaction  obtained 
Risk  Aversion-  concave  utility  function,  as  shown  below 
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Aside  on  Risk  Aversion/Tolerance 

(Jennings  &  Reichenstein,  2001) 

Pensions  considered  when  planning  retirement 
income . but  NOT  when  calculating  asset  allocation 

Pensions  and  investment  portfolio  generate  retirement 

funds;  why  not  consider  both  in  total  portfolio? 

Many  similarities  between  inflation-indexed  Treasury 
bonds  (TIPS)  and  military  retirement 

-  Linked  to  Consumer  Price  Index  (CPI) 

-  Backed  by  federal  government 

Suggest  treating  after-tax  present  value  as  a  “pseudo¬ 
bond”  in  total  portfolio 

Discounting  can  be  at  recent  TIPS  rates  (3%-5%)  or 
higher  personal  discount  rate  (18+%) 

Results  in  more  aggressive  (risk  tolerant)  portfolio  in 
active  investments  than  would  otherwise  result 
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NPV  of  Military  Retirement 


(Jennings  &  Reichenstein,  2001) 


Rank  at 
Retirement 

Years  of 
Service 

After-Tax 

NPV 

LTC 

20 

$726,674 

LTC 

22 

$802,690 

COL 

24 

$994,468 

COL 

26 

$1,096,490 

COL 

28 

$1,166,125 

COL 

30 

$1,205,255 

Assumptions: 

•Officer  currently  at  18  years  of  service 
•28%  tax  bracket 
•4%  TIPS  rate  /  inflation 
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“Pseudo-Bond”  Example 

(Nestler,  2007) 


40% 

Bonds 

Bonds 

60% 

f- 

Stocks 

Stocks 

Desired/Current 
Financial  Portfolio 
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Negative  Exponential  Utility 

U(w)  =  -e~cw,  c>  0 

•  Constant  Absolute  Risk  Aversion  (CARA)  -- 
no  “wealth  effect” 

•  Computational  tractability  advantage  over 
other  (log,  power)  utility  functions 

•  Analytical  solution  to  maximization  problem  is 
available  using  Certainty  Equivalent  (CE) 

•  CE  is  well-known  for  Normal  and  given  for 
VG-ICA  (Madan  and  Yen,  2004) 


Robert  H.  Smith 

SCHOOL  OF  BUSINESS 


Leaders  for  the  Digital  Economy 


©  2008  Scott  T.  Nestler 


Implied  Risk  Aversion  Coefficient 
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Portfolios  for  Comparison 


Model 

G  Fund 

F  Fund 

C  Fund 

S  Fund 

1  Fund 

VG-ICA 

(Daily) 

0% 

1% 

43% 

30% 

26% 

Riskless 

100% 

0% 

0% 

0% 

0% 

TSP 

“Market  Portfolio” 

39% 

6% 

39% 

8% 

8% 

L  2030 

16% 

9% 

38% 

16% 

21% 

L  2040 

5% 

10% 

42% 

18% 

25% 

k.  NOTE:  These  portfolios  are  created  with  returns  assumed  to  be 
Normally  distributed. 

Mr  Robert  H. Smith 
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Stochastic  Dominance 


•  Generalizes  utility  theory;  don’t  need  a  specific  utility 
function 

•  First-Order  Stochastic  Dominance  (FOSD) 

-  Assumes  only  monotonicity;  strongest  result 

-  A  FOSD  B  IFF  ^  (x)  >  ^  (x)?  VX 

•  Second-Order  Stochastic  Dominance  (SOSD) 

-  Also  assumes  risk  aversion 

-  A  SOSD  B  IFF  |  [FB(u)-FA(u)]du  >  0,  Vx 

—00 

•  Easy  to  test  with  empirical  data 

Robert  H.  Smith 
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Traditional  Risk  Measures 

•  Dispersion  Measures 

-Variance  (or Standard  Deviation) 

•  Treats  gains  and  losses  equally 

-  Semi-Variance 

•  Only  considers  observations  below  mean 

-  Mean  Absolute  Deviation  (MAD) 

•  Average  absolute  deviation  from  the  mean 

•  “Safety  Risk”  Measures 

-  Value-at-Risk  (VaR) 

-  Expected  Tail  Loss  (ETL) 

iSjt-  Robert  H.  Smith 
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Value-at-Risk  (VaR) 

-  “Expected  maximum  loss  over  a  fixed 
horizon  for  a  given  confidence  level” 

P(X  >  VaRx  (X))  =  /l 

-  Standard  risk  measure  for  past  12  years 

-  Does  not  reward  diversification 

-  Addresses  size  but  not  shape  of  tail 
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Coherent  Measures  of  Risk 

(Artzner,  Delbaen,  Eber,  &  Heath,  1999) 

•  Axioms  for  coherency: 

-  Translation  invariance  p(X  -  a)  =  p(X)  -  a 

-  Monotonicity  X  >  Y  =>  p(X)  >  p(Y) 

-  Sub-additivity  P(x  +  Y)  -  P(x)  +  P(Y) 

-  Positive  homogeneity  p(AX)  =  Xp(X) 

•  Variance:  not  monotonic  or  translation  invariant 

•  VaR:  not  sub-additive  in  non-Gaussian  world 

•  Other  measures  that  are  coherent  exist. 
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Conditional  VaR 


•  “Expected  value  of  all  losses  greater  than  VaR  for  a 
specified  X.” 

CVaRx  (X)  =  E[X  |  X  >  VaRx  (X)] 

•  Also  known  as  Expected  Shortfall  (Rockafellar  & 
Uryasev,  2001)  and  Tail  VaR  (Acerbi,  Nordio,  et  Al., 
2001) 

•  Accounts  for  size  and  shape  of  left  tail  but  ignores 
rest  of  distribution 


SCHOOL  OF  BUSINESS 


Leaders  for  the  Digital  Economy 


©  2008  Scott  T.  Nestler 


Classes  of  Weighted  VaR 

(Cherny,  2006;  Cherny  &  Madan,  2007) 

WVcRf(X )  =  f  CVaRx  (X)p(dx) 

"  LU?AJ 

•  Beta  VaR(a,(3) 

Hap  (clx)  =  B(J3  +  I,a-  /?)“'  xp  (1  -  xf^  dx,  x  e  [0, 1] 

-  Expectation  of  average  of  the  p  biggest  of  a 
independent  copies  of  portfolio  loss 

-  Faster  to  estimate  than  CVaR 

•  Alpha  VaR(a) 

-  Essentially  Beta  VaR  with  p=1 

-  Expectation  of  biggest  of  a  copies  of  portfolio  loss 

Robert  H.  Smith 
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mu(alpha,1) 


Effect  of  Alpha  and  Beta 


Effect  of  Alpha  With  Beta  =  1 


Effect  of  Beta  with  Alpha  =  50 


Can  allow  for  more  risk  by  decreasing  a  or  increasing  p 

Mr  Robert  H. Smith 
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Performance  (Reward-Risk)  Measures 


•  Sharpe  Ratio 

® £  w 
°x 

•  STARR  Ratio 

STARR  =  E ^ 

CVaR :x  (X) 

•  R-Ratio 

(Rachev) 

R  _  CVaR (-X) 
CVaRl2  (X) 

,f4fc  Robert  H. Smith 
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New  Portfolio  Performance  Measures 

(Nestler,  2007b) 

•  Similar  to  R-Ratio  but  use  Alpha-VaR 
and  Beta-VaR  in  place  of  CVaR 


•  AVaR-Ratio: 


AVR  = 


A  VaR gi  (-X) 
A  VaRai  (X) 


•  BVaR-Ratio:  g vr  =  BVaR«>^  X) 

BVaR^iX) 

Mr  Robert  H. Smith 

SCHOOL  OF  BUSINESS  © 2008 Scott T.  Nestler 


Leaders  for  the  Digital  Economy 


Monthly  Contribution 

•Assumes  saving  10%  of  base  pay  each  month  (median  for  TSP) 

$1,400 


$1,200 


1  21  41  61  81  101  121  141  161  181  201  221 

Month 


Robert  H.  Smith 

SCHOOL  OF  BUSINESS 


Leaders  for  the  Digital  Economy 


©  2008  Scott  T.  Nestler 


Realistic  Scenario:  Portfolio  Value 

(5000  sample  paths) 


Upside  Potential 

VG-ICA:  $  1,992,133 
L2040:  $  1,127,069 
L2030:  $  1,093,643 
TSP  MP:$  1,064,902 


Expected  Value: 

VG-ICA:  $418,381 
L2040:  $318,840 
L2030:  $313,936 
TSP  MP:$  310,247 
Riskless:$  259,642 


Mr  Robert  H. Smith  Day 
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Probability 


PDF  of  Discounted  Portfolio  Value 


-250000 

Robert  H.  Smith 
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Realistic  Scenario:  CDF  Comparison 
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NOTE: 

No  SD. 
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Realistic  Scenario:  Zoomed  CDF  Comparison 


95%  VaR 


tTy. 


i* 
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Realistic  Scenario 


&  Performance  Measures 


Risk  Measure 
(Jr  better) 

VG-ICA 

TSPMP 

L  2030 

L  2040 

Std  Dev 

$  168,885 

$  80,890 

$94,515 

$  105,525 

95%  VaR 

$43,382 

$44,146 

$  60,378 

$66,910 

95%  CVaR 

$  68,056 

$  54,789 

$  74,684 

$  82,783 

Alpha  VaR(50) 

$  77,575 

$  59,250 

$81,352 

$  87,754 

Beta  VaR(50,5) 

1  $  43,938 

$44,010 

$  60,203 

$  67,757 

Performance 
Measure  (t  better) 

VG-ICA 

TSPMP 

L  2030 

L  2040 

Sharpe  Ratio 

l  0.94  J 

0.63 

0.57 

0.54 

STARR  Ratio 

2.33 

2.90  | 

0.72 

0.69 

R-Ratio(.05,.05) 

8.87  | 

2.82 

4.10 

4.07 

AVR 

9.35  | 

5.41 

4.62 

4.67 

BVR 

[  1160  | 

5.05 

4.24 

4.12 

Mr  Robert  H. Smith 
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Traditional  Reward-Risk  Profile 


tTy. 


Risk  (Standard  Deviation) 


vs? 


Robert  H.  Smith 

jP  school  of  business 


r 
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Reward  (AVaR  (-X)) 


New  Reward-Risk  Profile 


Risk  (AVaR(X)) 

Mr  Robert  H. Smith 
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Some  Possible  Answers 

Why  might  I  be  more  risk  tolerant  than  I  currently  believe? 

-  Counting  military  or  government  pension  as  “pseudo-bonds”  could 
change  the  target  stock-bond  asset  mix. 

What  are  the  L  (Lifecycle)  funds?  How  are  they  constructed?  Why 
might  they  be  of  interest  (or  not)  to  me? 

-  “Set  it  and  forget  it”  funds  built  using  mean-variance  optimization  with 
returns  assumed  to  be  distributed  Normally. 

-  Depends  on  an  individual’s  level  of  interest  and  involvement. 

What  if  stock  and  index  fund  returns  are  not  normally  distributed,  as  is 
commonly  assumed? 

-  Possible  to  take  advantage  of  information  contained  in  higher  moments. 

How  does  choosing  reward-risk  measures  affect  optimal  TSP 
portfolios? 

-  Ability  to  capture  information  from  entire  distribution  is  useful. 

-  Need  to  do  further  work  on  optimizing  performance  measures  instead  of 
using  expected  utility. 

Robert  H.  Smith 

SCHOOL  OF  BUSINESS 


Leaders  for  the  Digital  Economy 


©  2008  Scott  T.  Nestler 


QUEST 
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